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Low-temperature thermodynamic properties of strongly interacting Fermi liquids with fermion 
condensate are investigated. We demonstrate that the spin susceptibility of these systems exhibits 
the Curie-Weiss law, and the entropy contains a temperature-independent term. The excessive 
entropy is released at the superconducting transition, enhancing the specific heat jump AC and 
rendering it proportional to the effective Curie constant. The theoretical results are favorably 
compared with the experimental data on the heavy fermion metal CeCoIns, as well as '^He films. 
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Theoretical understanding of strongly interacting 
Fermi systems, such as heavy-fermion metals, is challeng- 
ing P]. Conventionally, electrons in solids are classified 
as either itinerant or localized. The former have a Fermi 
surface, and their spin susceptibility x at low tempera- 
ture T follows the Pauli law x{T) — const, whereas the 
latter exhibit the Curie law x 1/T. The Curie law is 
observed in many heavy-fermion metals |^ Q |^ and 
is commonly attributed to the localized character of the 
/ electrons. However, in some of these materials, such 
as CeCoIn5 and PuCoGas [1,0, the low-temperature 
Curie law is immediately followed by a superconduct- 
ing transition, also associated with the / electrons 0- 
If the / electrons are localized, how can they supercon- 
duct? Actually, measurements of the Fermi surfaces of 
the heavy fermion metals by magnetic oscillations di- 
rectly demonstrate that the / electrons are itinerant, in 
agreement with band structure calculations @]. However, 
if the / electrons are itinerant, how can they exhibit the 
low-temperature Curie law? 

We show that these puzzles can be resolved within 
the Fermi-liquid theory if itinerant electrons form the so- 
called "fermion condensate" state 0, . The interplay 
between band-like and atomic-like behavior of electrons 
in solids is often treated on the basis of the Hubbard 
model and the dynamical mean- field theory |l j |. Heavy 
fermions are typically described by the Anderson-Kondo 
lattice models of coupled itinerant and localized electrons 
originating from different orbitals, sometimes using the 
two- fluid description [l^ . However, given the experimen- 
tal evidence that the / electrons are itinerant in some 
heavy-fermions metals, here we study a conceptually sim- 
pler model where all electrons are itinerant. Our goal is 
not to present a detailed, material-specific description, 
but to illustrate general ideas also applicable to other 
puzzling Fermi systems, such as '^He films 

Let us consider a system of itinerant electron quasipar- 
ticles characterized by dispersion Sp, where e is energy 
measured from the chemical potential, and p is momen- 



tum. The spin susceptibility xo per one electron is 



Xo 



dn{£p) 



n{p) [1 - n{p)] dvp 



(1) 

where dvp = 2d^p/ p {2TThY is the volume element in 3D 
momentum space, p is the electron concentration, and 
/ie is the magnetic moment of electron in solid. In a 
simple case, p^. is equal to the Bohr magneton pB^ but, 
for heavy fermions, Pe niay also contain a contribution 
from the orbital angular momentum, as discussed later 
in the paper. The occupation numbers n(p) are given by 
the Fermi distribution function 



n{p) = [l + exp(ep/r)]- 



(2) 



where the Boltzmann constant ks is set to 1. In ordinary 
Fermi liquids, Ep « vp{p — pp), where pp is the Fermi 
momentum, and vp is the Fermi velocity. In this case, 
the integral in Eq. ^ is proportional to T, and XoiT) = 
plNf)/p = const, where Nq = p'^/n'^vp is the density of 
states at the Fermi level. Accounting for the spin-spin 
interaction amplitude go modifies Eq. (QJ via the Stoner 
factor: x = Xo/(l - 9oXo)- 

The quasiparticle dispersion Ep is affected by the Lan- 
dau interaction function fL{p,p') — 5ep/5n{p'). In gen- 
eral, fL{p,p') is a functional of the occupation num- 
bers n{p), but here, for the sake of illustration, we take 
Il{PtP') as a given function. Then, Sp is related to the 
bare dispersion as 



fL{p,p')n{p')dv. 



p' ■ 



(3) 



The dispersion Ep and the occupation numbers n{p) are 
obtained by solving Eqs. Q and sclf-consistently. 
When the interaction is weak, Eq. ||2J) merely renor- 
malizes the Fermi velocity. However, when the interac- 
tion strength exceeds a critical value, the minimum of 
the total energy at T = may be achieved in a radi- 
cally different state with the fermion-condensate 0- In 
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this state, the quasiparticle spectrum is completely flat 
Ep = at T = in some region of momentum space, 
where the occupation function n* (p) continuously inter- 
polates between and 1 |^. The increase of kinetic 
energy in this state is compensated by the decrease of 
interaction energy. Nozieres 0| demonstrated that, at 
low T 7^ 0, the momentum occupation function remains 
the same n{p) — n^,{p) in the domain occupied by the 
fermion condensate. Thus, the self-consistent dispersion 
Ep becomes temperature-dependent 



£p = T In 



1 - "*(P) 



(4) 



as follows from the inversion of Eq. 10). The group ve- 
locity dsp/dp in Eq. is proportional to T, which 
generates a sharp peak in the density of states with the 
height proportional to l/T and results in unusual ther- 
modynamic properties discussed below. Measurements 
of magnetic oscillations in the heavy-fermion metals in- 
deed show enormous flattening of £p relative to the band- 
structure calculations ,8j. 

This qualitative analysis has been confirmed by ana- 
lytical and numerical solutions of Eqs. Q and lO for 
various interaction functions f^. As an example, in Fig. 
n we show n{p) and Sp numerically calculated for three 
different temperatures for a toy model with an isotropic 
parabolic dispersion, characterized by the bare mass m 
and the bare Fermi energy — p^/2m. The interac- 
tion function fL{q), where q = p — p', was chosen to 
be fL{q) = A/{[1 - (g/2pi.)2]2 + with P = 0.48 and 
A = Sppm/ir'^. Panel (a) shows that n{p) markedly dif- 
fers from a step function and does not depend on tem- 
perature in the interval [pi,p2]- Panel (b) demonstrates 
that Sp changes with T, but the inset shows that the 
ratio Sp/T is T-independent in the interval [pi,p2], m 
agreement with Eq. (@J. Fig. illustrates that a self- 
consistent solution of quite conventional equations (0) 
and for a rather generic, non-singular interaction 
function does generate the fermion condensate. Simi- 
lar results were found for other isotropic and crystal lat- 
tice models, where p is quasimomentum in the Brillouin 
zone: see Ref. jl^l and references therein. The results 
are robust and do not depend significantly on model de- 
tails. Although the toy model utilized for the calculations 
shown in Fig. ^ is not necessarily realistic for specific 
materials, the fermion condensate formation is a generic 
process, and only the fermion condensate parameters, not 
the model details, matter for observable quantities. 

Now let us discuss observable manifestations of the 
fermion-condensate state in detail. Substituting the 
temperature- independent occupation function n^{p) of 
the fermion condensate into Eq. Q i we find a Curie con- 
tribution to the spin susceptibility |16l] 



even though all electrons are itinerant. The effective 
Curie constant in Eq. (01 is reduced by the dimension- 
less parameter k relative to the Curie law xo = fJ'l/T of a 
non-degenerate Fermi gas at high temperatures T > e'p. 
The second term x in Eq. (O comes from integration 
outside of the fermion-condensate domain [pi , P2] in Eq. 

This term is less singular than the Curie term, 
which dominates at low T. Accounting for the spin in- 
teraction amplitude go generates the Curie- Weiss law 
x(T) « ^1n/{T — Qw) with the Weiss temperature 
©VK = So^Me- The numerically calculated Xq{T) for the 
same model as in Fig. ^ is shown by the solid line in 
Fig. [2^.. The finite value of the product xoT in the limit 
T — > indicates the Curie behavior at low temperatures 
and gives the value k « 0.1 in this model. In a wide tem- 
perature range, Xq{T) shown in Fig. |2K does not strictly 
follow the Curie law because of x{T)- 

In ■^He films, the low-T Curie constant is about 4 
time lower than the high-T one, as shown in Fig. 1 of 
Ref. 01, which gives k ~ 0.25 in this case. To eval- 
uate K for the heavy fermion metals from Eq. (jSJ, we 
need to know the magnetic moment of / electrons, 
which has spin and orbital contributions. In a free atom, 
/X — gLfJ-sJ^ where is the Lande factor, and J is 
the total angular momentum. The crystal field lifts de- 
generacy between energy levels with different projections 
Jz and causes magnetic anisotropy. In CeCoIus, there 
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n^{p)[l~n^{p)]dvp , (5) 
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FIG. 1: Quasiparticle properties of a system with the fermion 
condensate, plotted vs. p/pF for three temperatures T: (a) the 
occupation numbers n(p), (b) the single-particle spectrum Sp 
in the units of e''p — p%/2m, (inset) the ratio Sp/T. 
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is one /-electron with J = 5/2, and the c axis is the 
easy magnetic axis. Thus, the lowest energy levels have 
Jz — ±5/2, and the effective magnetic moment along the 
c axis is /ig = ^-bQlJz — where = 5/2 and 

= 0.857 for L = 3, S* = 1/2, and J = 5/2. As shown 
in Fig. 3 of Ref. 0], the low-T Curie law in CeCoIns is 
the most pronounced for the easy axis c with the Curie 
constant 0.2/i^. This value is much smaller than /ig, and 
we find from Eq. © that k = 0.2^% / {2. U^b)'^ = 0.044. 

The entropy S per one particle for an ensemble of 
fermion quasiparticles is given by the formula 



{n{p)\iin{p) + [l-n{p)]lTi[l-n{p)]}dvp . (6) {p-P2f for Ep > 0, which gives C cx T^/^ and x oc T 



as pressure P, so that dS^/dP ^ 0. 

In ordinary Fermi liquids, the specific heat C — 
TidSjdT) is proportional to T. The straight line of 
slope 1/2 in Fig. |3i shows C/tt^ and XoT/2,ijlI vs. T/e^, 
for a Fermi gas. The fermion condensate does not con- 
tribute to the specific heat, because its entropy 5* is T- 
independent. So, the main contribution to C comes from 
the regions in momentum space where Sp interpolates 
between the fermion condensate and the regular disper- 
sion (see Fig. When the fermion-condensate domain 
is small, Ep has inflection points [l^ Il7| , and the lead- 
ing term in Ep is proportional to {p — piY for £p < and 

-2/3 



In ordinary Fermi liquids, the integrand in Eq. ^ differs 
from zero only in a narrow vicinity of the Fermi sur- 
face, so S* = Ttt'^ /vppp, and S ^ when T ^ 0. In 
contrast, in the fermion-condensate state with the occu- 
pation function n^{p), the integrand is nonzero in a finite 
region, so the entropy has a temperature-independent 
contribution 5*. Extrapolating the normal-state S{T) 
from the inset in Fig. 2 of Ref. Q to T ^ 0, we esti- 
mate that 5", « 0.1 In 2 in CeCoIn5. We see that 5* is 
finite, but much smaller than In 2 expected for an en- 
semble of localized spins 1/2. Because it is proportional 
to the momentum-space volume occupied by the fermion 
condensate, 5** may depend on external parameters, such 
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FIG. 2: Thermodynamic properties of a fermion-condensate 
system. (a) The spin susceptibihty xoT/3fie (solid line) 
and the specific heat C/tt^ (dashed line) vs. temperature T. 
The inset: the same variables in the log-log scale, (b) The 
Sommerf eld- Wilson ratio (solid line) vs. T. The dotted lines 
in panels (a) and (b) correspond to the Fermi gas. 



The dashed line in Fig. |2K shows the numerically calcu- 
lated C{T) for the same model as in Fig.^ which indeed 
exhibits a sublinear power law. Because the calculated 
XqT and C have different temperature dependences, the 
Sommerfeld- Wilson (SW) ratio R^^w = ^^'^XoT/SCpl in- 
creases with decreasing T, as shown by the solid line in 
Fig. and observed in heavy fermions . This is in 
contrast to ordinary Fermi liquids, where R^g^ = 1, as 
shown by the horizontal line in Fig. |2b- Notice that the 
Stoner factor is not included in our definition of Pgyj/. 

Although the excessive entropy S*, of the fermion 
condensate does not contribute to the specific heat, it 
produces an enormous enhancement of the thermal ex- 
pansion coefficient a = dV/dT = —dS/dP and the 
Griineisen ratio T — a/C [l9|. In ordinary Fermi liq- 
uids, S (X T, thus a oz T vanishes at T — > 0, and 
r(T) = const. In contrast, for the fermion condensate, 
the derivative dS^/dP is T-independent, so a has a finite 
value at T — > 0. Experiment [2(J shows that a is indeed 
temperature-independent at low T and exceeds typical 
values for ordinary metals by the factor of 10^-10'''. With 
a — > const and C(T) — > 0, the Griineisen ratio F = a/C 
diverges at low T, which is observed experimentally |2l| . 

However, the existence of the residual entropy S^, at 
T — > contradicts the third law of thermodynamics (the 
Nernst theorem). To ensure that 5* = at T = 0, lo- 
calized spins order magnetically due to spin-spin inter- 
action. Similarly, a system with the fermion condensate 
must experience some sort of a low-temperature phase 
transition eliminating the excessive entropy 5*. Here we 
focus on the second-order phase transition to a supercon- 
ducting state . The progressive increase of the fermion- 
condensate density of states with decreasing temperature 
facilitates superconducting instability in one of the pair- 
ing channels: s, p, d, etc. Elementary excitations in a 
superconductor are the Bogolyubov quasiparticles, whose 

spectrum Ep = £p + has the energy gap Ap. The 
entropy of a superconductor is given by Eq. © with 
^(£p) ^ /(-^p)j where f(Ep) are the occupation num- 
bers of the Bogolyubov quasiparticles. Because of the 
energy gap, / ^ 0, and so 5 ^ at T ^ 0, thus the 
Nernst theorem is satisfied. However, in order to release 
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the excessive entropy S'*, the specific heat jump AC at 
the transition temperature Tc is enhanced. 

The specific heat Cg of a superconductor is Cg = 
J dvp Ep df {Ep) / dT . Taking the difference between Cs 
and Cm the specific heat in the normal state, we find the 
specific heat jump at Tc- 

AC = C, - C„ = 1 ^n{p) [1 - n{p)] dv^ . (7) 

In the BCS theory, dA^/2TcdT = -47rV7C(3) « -4.7, 
where C is the zeta function. (For d-wave pairing, the 
number is different, but the results do not change signif- 
icantly.) Comparing Eqs. ^ and {Tj), we find 

AC = 4.7TeXo/Me = 1-43 C„4h' • («) 

Eq. © shows that the specific heat jump AC can be 
expressed in terms of either Xo(7c) or Cn{Tc). For or- 
dinary Fermi liquids, where the Sommerfeld- Wilson ra- 
tio bI'sw ~ 1' © reproduces the familiar BCS rela- 
tion AC/C„ ~ 1.43. However, in the fcrmion condensate 
state, changes with temperature, so AC /Cn does 

not have a universal value. Since xo ^ ^/T, the specific 
heat jump ijSJl is not proportional to Tc, but is related to 
the fermion condensate parameter k in the Curie law (O 

AC « 4.7 k. (9) 

Thus, the ratio AC/C„ can be very high when Tc is low, 
because C„ ^ at T — * while AC is finite. 

Let us apply this analysis to CeCoIns, where Tc = 2.3 
K. In this material AC/C„ ~ 4.5 is substantially 
higher than the BCS value, in agreement with our ar- 
guments. Using the value k « 0.044 evaluated from the 
Curie law, we estimate the r.h.s. of Eq. ^ as 0.21. This 
is about a half of the experimental value of AC w 0.42 
per electron measured in Ref. |2^ . However, as shown in 
Fig. 3 of Ref. [H , the Curie term constitutes only about 
a half of the spin susceptibility at Tc- Thus, using the 
total susceptibility, we find that the relation © between 
AC and xo is satisfied. In PuCoGas, we also attribute 
the Curie law, followed by a superconducting transition 
at Tc = 18.5 K j^, to the fermion condensate. How- 
ever, quantitative estimate of k is difficult in this case, 
because plutonium has five / electrons. Interestingly, if 
plutonium is replaced by uranium, the resulting material 
UCoGas does not exhibit the Curie law and does not have 
superconducting transition [^. This is surprising from a 
conventional point of view, where "itinerant" electrons in 
UCoGas should be more susceptible to superconducting 
pairing than "localized" electrons in PuCoGas. However, 
if the fermion condensate does not form in UCoGas, so 
that there is no Curie law, then the density of states is 
not enhanced, and superconductivity is not facilitated. 

In conclusion, we have shown that strongly interact- 
ing Fermi liquids can form a fermion-condensate state. 



where quasiparticle dispersion Sp is flat at the Fermi 
level. Their magnetic susceptibility Xo{T) exhibits the 
Curie- Weiss law with the effective Curie constant reduced 
by the fermion-condensate parameter k. The entropy 
has the temperature-independent term S"* (estimated as 
5* « 0.1 In 2 per electron in CeCoIns), which greatly in- 
creases the thermal expansion coefficient a — —dS/dP 
at low T. The excessive entropy S"* is released below the 
superconducting transition temperature Tc, which dra- 
matically reduces a and enhances the specific heat jump 
AC/C„, as observed in CeCoIns 0, The univer- 

sal relation ((SJ between AC and TcXo{Tc) can be tested 
experimentally by checking whether the both quantities 
change proportionally upon variation of externalparam- 
eters, such as pressure or chemical substitution |7ll24l|. 
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